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Context. The Solar radiation field may break apart ultra high energy cosmic nuclei, after which both remnants will be deflected in the 
interplanetary magnetic field in different ways. This process is known as the Gerasimova-Zatsepin effect after its discoverers. 
Aims. We investigate the possibility of using the detection of the separated air showers produced by a pair of remnant particles as a 
way to identify the species of the original cosmic ray primary directly. Event rates for current and proposed detectors are estimated, 
and requirements are defined for ideal detectors of this phenomenon. 

Methods. Detailed computational models of the disintegration and deflection processes for a wide range of cosmic ray primaries in 
the energy range of 10" to 10^" eV are combined with sophisticated detector models to calculate realistic detection rates. 
Results. The fraction of Gerasimova-Zatsepin events is found to be of the order of 10"^ of the cosmic ray flux, implying an intrinsic 
event rate of around 0.07 km"^ sr"' yr~' in the energy range defined. Event rates in any real experiment, however, existing or under 
construction, will probably not exceed 10"^ yr"' . 
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1. Introduction 

The mass composition of very high energy cosmic rays provides 
important information on their acceleration mechanisms and the 
compositions of their sources. Usually, however, it is only pos- 
sible to make statistical, model-dependent estimates of the pri- 
mary particle types from an ensemble of showers. Primary com- 
positions are then derived from the abundances of species com- 
ponents in the air showers considered ( Antoni et al.|2002 l. 

An alternative mass determination makes use of the Gera- 



simova-Zatsepin effect (Zatsepin 1951 ; Gerasimova & Zatsepin 
|1960| l. In this scenario, one relies on the fact that atomic nuclei 
have a chance of undergoing photodisintegration in the Lorentz 
boosted Solar radiation field before arriving at Earth, splitting 
the nucleus into two parts. Due to the different rigidities of these 
two fragments, the deflection in the interplanetary magnetic field 
will be different, resulting in two separate air showers at some 
spatial separation, but arriving essentially at the same time and 
from the same direction ( Medina-Tanco & Watson|1999 Epele 
|et al.|[T999| ). To our knowledge, no experimental detection of a 
Gerasimova-Zatsepin event has ever been reported. 



2. The Gerasimova-Zatsepin process 

The photodisintegration probability rjz for a nucleus approach- 
ing the Solar system to undergo photodisintegration has been 
investigated tho r oughly by |Zatsepin| (|1951|l, [Gerasimova & 
Zatsepm (T960|, [Medina-Tanco & Watsonr ( |1999[ ) and [Epele 
et al. (1999). It can be calculated by integrating its path 



length against photodisintegration over its trajectory. The pho- 
ton energy as seen from the cosmic ray's comoving frame is 



Lorentz boosted by a factor 2'>'cos^(q'/2), where y is the cos- 
mic ray's Lorentz factor and a is the angle between the prop- 
agation directions of photon and particle in the heliocentric 
frame ( Gerasimova & Zatsepin 1960). Different photodisintegra- 
tion reactions are possible, but by far the most likely reaction to 
occur, is those in which one proton or neutron is knocked out of 
the nucleus ( Karakula & Tkaczyk|1993 1. 

After disintegration, the charged remnants will be deflected 
in the interplanetary magnetic field. Since the mass/charge ratio 
will generally be different for the two fragments of the disinte- 
grated nucleus, so will the amount of deflection be. The shape 
and strength of the magnetic field surrounding the Sun is quite 
complicated. Akasofu et al. (1980i have constructed a three- 
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dimensional model which consists of four components: (i) the 
Solar dipole magnetic field; (ii) a large number of smaller mag- 
netic dipoles located along an equatorial circle just inside the 
Sun; (iii) the field of the poloidal current system generated by 
the Solar unipolar induction; and (iv) the field of an extensive 
cuiTent disc around the Sun, lying in the ecliptic plane. The in- 
fluence on the deflection of cosmic particles is dominated by the 
latter two components, as their contributions is larger at greater 
distances. 

Given the discreteness of the masses of the remnants and the 
linear proportionality between a remnant's mass and its energy 
(assuming single-nucleon emission), the mass number A of the 
original disintegrated particle can simply be determined by esti- 
mating the energies of the primaries of the two showers ( jEpele| 
|etal.|1999D : 

A.^l±^, (1) 

£i ' 

where £1 is the energy of the less energetic shower. 
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3. Detection of Gerasimova-Zatsepin events 

Identifying a Gerasimova-Zatsepin pair as such requires both 
showers to be seen by a cosmic ray detector. Cyclotron radii 
for cosmic rays at energies above 10'*' eV and magnetic field 
strengths < 10"^ T are very large compared to the size of the 
Solar system, allowing us to take both remnants' arrival direc- 
tions to be equal to each other and to the original arrival di- 
rection. In order to calculate the Gerasimova-Zatsepin detection 
aperture for a given cosmic ray detector, let us define the sepa- 
ration resulting from different amounts of deflection of the two 
showers as the vector d = (6\\,6j_) between the two tracks upon 
impact, transverse to the arrival direction ((po, Oq) in the Solar ref- 
erence frame. Let 6\\ be the component in the ecliptic plane and 
((po, 6o) - (0, 0) in the direction of the Sun. 

For an accurate description of a detector's aperture, it is nec- 
essary to incorporate the angle at which the detector is hit by the 
cosmic ray particles. The angles due to daily and yearly phase 
of Earth, ^d and ^y respectively, together with the latitude b and 
longitude / of the detector, fix the orientation of the detector with 
respect to the sky. The distance between the two showers 6' de- 
pends on these four angles, and affects detection rates in a non- 
trivial way. As the ratio 6' /6 may easily exceed a factor of 2, it is 
clear that projection effects cannot be neglected in our analysis. 

Whether both air showers are actually detected, depends 
on the detector geometry as laid out on Earth. Let us define a 
detector-specific function ^(6) which describes the probability 
of detecting the second shower event for a given separation vec- 
tor 6, under the assumption that the first shower is detected. The 
effective aperture A can now be calculated by integrating over ^ 
over the course of a year. We are also taking into account the de- 
tector's angular sensitivity w as a function of the zenith angle 6 
as observed by the detector: 
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Fig. 1. Absolute Gerasimova-Zatsepin energy spectrum (solid 
line), and upper limit for Earth-based detectors, i.e. events with 
\6\ < 2Rg) (dashed line). Also shown is the fraction of the overall 
galactic hadronic cosmic ray flux Jq in the bottom panel, show- 
ing that 10"^ < T/Gz < 10""^ in this energy range. 



tor for particles with energies greater than E is given by 
Ogz(£)= r 2/z(£') fjiziE',cf>o,eo)A(6,E',<Po,0o) 
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where 5 o is the total area covered by the detector, and < w < 1 . 
The factor 1 /n serves to normalise to all sky visibility. 

The absolute particle fluxes for various primary nuclei are 



where /dc is the duty cycle of the detector, which is a constant 
factor in case of surface scintillators, but may depend on 0o 
and 00 for example for air fluorescence detectors, as they can- 
not observe during the day. 



estimated from the model presented by Horandel (2003 i, which 4. Results 



Jz(E) = 7o,z 



E 



yz 



EpZ 



(3) 



Jz(E) are the contributions of a species Z to the cosmic ray 
spectrum, 7oz and yz are constant factors for each species, 
Eo = 10'^ eV, Ep = 4.5 ■ 10'^ eV, ji = 1.9 and 72 = l-l- For 
E > 10'*^ eV, this model underestimates the number of cosmic 
ray particles of very low mass in the spectrum. However, pho- 
todisintegration cross sections at the these energies are too low 
to be of consequence for very light nuclei. The total hadronic 
cosmic ray flux is 

" (4) 



y(£) = ^7z(£), 



where the summation runs over all cosmic ray particle species. 
In our case 2 < Z < 92, as protons will not contribute to the 
Gerasimova-Zatsepin flux at all. 

If TjziE, (pi), 0q) is the probability for a nucleus of species Z 
and energy E to undergo photodisintegration along its trajectory, 
then the final Gerasimova-Zatsepin event rate for a given detec- 



To calculate realistic values for rjz, a numerical model was 
constructed. Calculations were carried out for primary cosmic 
ray species from "^He to ^■'^U, with energies ranging from 10'^ 
to 10^° eV, the region where the photodisintegration cross sec- 
tion is highest. The obtained average values of rjz for Fe, O 
and He are in line with earlier findings ( jEpele et al.|| 19991 
Medina- Tanco & Watson|1999 l. The partial contribution of the 
heavier nuclei to the Gerasimova-Zatsepin spectrum is larger 
than one might expect, due to their high overall value of rjz- 

By multiplying each species' disintegration probability by 
its partial flux according to Eq.13] the total intrinsic Gerasimova- 
Zatsepin flux Jgz{E) - 2z JziE) J T]z cos 60 d0Q dcpo is obtained. 
Fig. [1] shows this flux as a function of energy. The solid line 
represents the absolute total flux by counting all disintegration 
events. For reference, the flux relative to the integral cosmic ray 
spectrum J{){E) is also drawn in the bottom panel, showing a 
maximum disintegration probability of ?7gz - 10"** near E ^ 
1.5 ■ 10'^ eV. The dashed line was obtained by disregarding any 
event with a separation larger than one Earth diameter. This line 
sets a hard upper flux limit for any Earth-based detector. Notice 
that events with these very high separations primarily occur in 
the lower energy end of the spectrum: this makes sense, as the 
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Fig. 2. Probability distribution of separations for species of "^He, 
'^O, ^'^Si, ""^Fe and ^""^U. Thick lines are for proton emission, thin 
lines denote neutron emission. Shown are expected separations 
for a primary of £ = 10"* eV; for other energies, multiply 6 
by 10'** eV/£. 



separation of a disintegrated cosmic ray pair is expected to be 
proportional to the inverse of its energy. 

The disintegration probability strongly depends on the ar- 
rival direction. The disintegration process favours arrival direc- 
tions close to the Sun, as higher integrated photon densities boost 
the number of disintegrations over trajectories from this direc- 
tion. 

Given the complexity of the magnetic field in the Solar sys- 
tem, a numerical equivalent of the field was implemented and 
disintegrated particles were propagated accordingly. Particle tra- 
jectory deviations are largest for directions near the Sun: too 
large, in fact, to be detected. Therefore, counterintuitively, high- 
est event rates for any realistic cosmic particle observatory are to 
be expected on the night side of the sky, even though fluxes are 
lower there. 

To illustrate the dependence of the separation distance on the 
various parameters, the distribution of 6 is given in Fig. l2]for se- 
lected species. All curves are normalised so that their respective 
logarithmic integrals equal unity; separate lines are drawn for 
disintegration reactions involving proton and neutron ejection. 
The separations shown are for a primary of E - 10'** eV. Apart 
from statistical deviations, all curves are identical when separa- 
tions are shifted to the left by a factor A|Zi /Ai - Z2/A2I, peaking 
at 40 km. Since 6 oc £"', multiplying the separation value by 
10"* eV/E yields the correct separation at other E. We may now 
parameterize the expectation value for the separation as 



{6} = 4A 






A2 



10''' eV 



km. 



(6) 



The overall remaining shape of the separation distribution is a 
result of the magnetic field shape and strength and the disinte- 
gration distance R from Earth only. 

Let us now estimate event rates for selected air shower exper- 
iments: the existing Pierre Auger Observatory (Abraham et al. 
|2004| l and the Lofar radio telescope (Falcke et al. 2006|l, which is 
under construction. Auger and Lofar are both surface detectors, 
sensitive to secondary shower particles (Auger) and radio signals 
produced in the shower (Lofar). Auger is a dense array, covering 
a continuous area with a total aperture of A = 4.7 ■ lO-' km^ sr 
Lofar's geometry is sparser, consisting of interconnected smaller 
stations with no detectors in between. Though Lofar's aperture 
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Fig. 3. Expected integrated Gerasimova-Zatsepin detection rates 
for Pierre Auger (solid lines) and Lofar (dashed). Four lines are 
drawn for each detector, representing, from top to bottom, the 
theoretical upper limit for a detector of the given area; applying 
separation cut-off's regarding detector geometry; applying loose 
energy cuts; and applying strict energy cuts (see text for further 
explanation). 

is much smaller at A = 2.2 ■ 10^ km^ sr, it is able to recon- 
struct showers with a much lower energy. Both detectors were 
modelled numerically to obtain accurate values for the effective 
aperture and Gerasimova-Zatsepin event rate. For each detector, 
simulations were carried out to make predictions for the final 
event rates <^cz(E) according to four scenarios: 

1 . As a simple first step, we can set a hard upper limit by taking 
the observatory's total aperture as effective cross section, im- 
plying Ocz(E) - 0.45A km"^ sr"' yr"' for energies between 
10'^ and 10^'^ eV. This approach means that every event has 
nonzero probability of being detected, regardless of its rem- 
nants' separation. 

2. A more realistic estimate is obtained by applying the aper- 
ture function A according to Eq. |2] In this way, we include 
projection effects as a result of the detector's orientation. We 
also apply a lower limit 6mm to the separation distance; this 
is the minimum separation at which the detector can disen- 
tangle two showers. 

3. In this scenario, a further restriction is applied by imposing 
a lower limit E^m on the more energetic shower For the less 
energetic shower, an energy down to a tenth of this limit is al- 
lowed. This approach is justified by the possible implemen- 
tation of a triggering system in which data for the entire de- 
tector array is stored for each trigger, allowing one to check 
for coincidences at a later time. 

4. By applying a strict energy cut, demanding that both show- 
ers exceed the threshold energy, a less sophisticated trigger 
suffices. This scenario is a pessimistic assumption for detec- 
tors which are not optimised for Gerasimova-Zatsepin pair 
detection. 

Both detectors' ^(6) functions were generated numerically 
from their geometries. The values for E^m - 10'^ eV for Auger 
and 10'^ eV for Lofar were interpreted as Gaussian error func- 
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Fig. 4. Required aperture for detection of a single Gerasimova- 
Zatsepin event per year as a function of £„!„ for a circular detec- 
tor at 45° latitude, /dc = 1, ^min = 60° and Smm = 1 km. 



tions. Maximum zenith angles are 60° in the case of Auger and 
80° for Lofar. 

Derived event rates for each scenario are presented in Fig. [3] 
This figure shows separate lines, from top to bottom, for each 
scenario. When the observatory's lower energy limit is taken into 
account according to scenario 4, event rates plummet to levels 
between 10"^ and 10""^ yr"'. This figure effectively dismisses 
any possibility of successful Gerasimova-Zatsepin pair record- 
ing. Scenario 3 is certainly a possibility in the case of Auger, as 
information for all surface detector tanks is stored when a big 
event is seen. Even in scenario 3, however, event rates are ex- 
pected to stay below 10"^ yr"'. 

Lofar's detection area is about 36 times smaller than that of 
Auger, and consequently it has a much lower intrinsic Gerasi- 
mova-Zatsepin flux, not exceeding 0.3 yr"' . Lofar's energy limit 
is 10 times lower, however: this easily compensates the lack of 
area, as the slope of the Gerasimova-Zatsepin event rate is ap- 
proximately oc E^^ at higher energies (see Fig.Bj). Still, expected 
event rates do not exceed 5 ■ lO"'* yr"' in scenario 4. Comparing 
event rates from simulations of the full detector and the central 
core only, shows that no significant contribution to pair detec- 
tion is to be expected from the relatively small outer stations. 
The increased rate in scenario 3, probably the maximum achiev- 
able rate for Lofar, is still quite insignificant at 2 ■ 10"^ yr"'. 
The chances of the corresponding trigger algorithm being im- 
plemented are slim at best, because the Lofar telescope will not 
be a dedicated cosmic ray detector, but an experiment shared 
with other astronomical observations. In practice, data band- 
width limitations would probably not allow a trigger to be com- 
municated to every antenna in the array, except for very energetic 
events which produce a negligible Gerasimova-Zatsepin rate. 

Because the expected event rate for a cosmic ray detector 
is approximately proportional to SqE^^, it pays off to invest in 
lowering a detector's energy threshold instead of focusing on 
collecting area if maximising Gerasimova-Zatsepin event rates 
is intended. Still, any detector would have to be huge to receive 
one single detectable event per year or more. Estimated required 
aperture sizes are shown in Fig. 4]for a detector with a Gaussian 
lower energy limit £■,„!„ with full duty cycle and 60° zenith limit. 
At Auger's energy limit, the required aperture in scenario 3 
would be ~ 90 times that of Auger; at Lofar's energy limit, it 
would be ~ 50 times that of Lofar Even at low threshold en- 
ergies, where an aperture of only ~ 1.4-10^ km^ sr would be 



required, the cost of constructing and operating such an array 
would easily amount to several billions of euros. 

For surface detectors, there might be other points to consider. 
Even if a cosmic ray observatory was hit by a detectable Gera- 
simova-Zatsepin pair, it is very much the question whether such 
an event would meet the quality criteria. The occurrence of two 
overlapping or nearly overlapping showers within a very short 
time window might prevent proper reconstruction, and the event 
would be discarded as noise. 

Fluorescence detectors will not produce more favourable Ge- 
rasimova-Zatsepin detection rates, as their duty cycles are typ- 
ically only 0.1, increasing the aperture needed by an order of 
magnitude over surface detectors. Moreover, the fluorescence 
detection technique does not work well for low energies, where 
rates are highest: typically, a 10'^ eV air shower's fluorescence 
signal can be seen no more than 10 km away from a fluorescence 
telescope. This also discards satellite missions such as euso as 
possible detectors. 

On a side note, an interesting deviation of the pure cosmic 
ray spectrum as seen by current experiments can be observed. 
We have shown that a fraction of ~ 10"^ of all cosmic particles 
that aiTives at Eaith does so as a Gerasimova-Zatsepin pair. In 
practice, these events go unnoticed as such and are registered as 
normal cosmic particles. This means that, at any given energy, 
a small fraction of the events is to be attributed to a primary 
of much higher energy. The actual spectrum then has a spectral 
index of the order of 10"^ less steep than currently assumed; this 
factor is too small to perceive given current data eiTor margins. 

5. Conclusion 

We have used a set of simulations to calculate in detail the rate 
of very high energy cosmic ray particles breaking up in the Solar 
magnetic field and the probability distributions of the separations 
of their remnants. Additionally, we have used accurate detector 
models to estimate realistic detection rates for this phenomenon. 
We have shown that current experimental setups, including the 
PieiTe Auger Observatory, by far lack either the energy sensitiv- 
ity or the area to produce any significant amount of detections 
of the kind, and would detect only a fraction of the Gerasimova- 
Zatsepin flux predicted by Epele et al. ( 1999 1 and Medina- Tanco 
& Watson (1999 ). Consequently, the prospects for any future ex- 
peiiment detecting the Gerasimova-Zatsepin effect are negligi- 
ble. 
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